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The locations of the peaks of the CMB spectrum are sensitive indicators of cosmologi- 
cal parameters, yet there is no known analytic formula which accurately describes their 
dependence on them. We parametrize the location of the peaks as Im = lA{m — fm), 
where Ia is the analytically calculable acoustic scale and m labels the peak number. 
Fitting formulae for the phase shifts ipm for the first three peaks and the first trough 
are given. It is shown that in a wide range of parameter space, the acoustic scale Ia 
can be retrieved from actual CMB measurements of the first three peaks within one 
percent accuracy. This can be used to speed up likelihood analysis. We describe how 
the peak shifts can be used to distinguish between different models of dark energy. 



1 Introduction 

The locations of the peaks and troughs of the CMB 
anisotropy spectrum can serve as a sensitive probe of cos- 
mological parameters [1, 2, 3, 4, 5, 6]. 

There are however many processes which contribute to 
the final anisotropies, and these must be calculated from 
systems of coupled partial differential equations [7]. As 
such it is not possible a priori to derive an accurate analytic 
formula for the peak locations. There exists a numerically- 
obtained estimate of the location of the first peak [8] for 
a universe with no cosmological constant, namely l\ ~ 

1/2 

200 VLm ' . This was recently extended to universes with 
A 7^ 0, by perturbing around the A = value [9], but 
holding all other parameters fixed. In this work, we calcu- 
late the locations of the first three peaks as a function of 
several cosmological parameters, including universes with a 
large dark energy component. We show how these results 
can be used to extract cosmological information about, for 
instance the history of quintessence, from just a handful 
of CMB data points and also to speed up multi-parameter 
likelihood analysis. 

Before last scattering, the photons and baryons are 
tightly bound by Compton scattering and behave as a fluid. 
The oscillations of this fluid, occuring as a result of the bal- 
ance between the gravitational interactions and the pho- 
ton pressure, lead to the familiar spectrum of peaks and 
troughs in the averaged temperature anisotropy spectrum 







l\ (estim.) 


li (numeric.) 


% error 


0.4 


0.6 


296 


219 


35 


1.0 


0.0 


269 


205 


31 



Table 1: Values of the location of the flrst peak h estimated 
by h ~ I A and calculated numerically via CMBFAST [7]. 
The intuitive model clearly does not describe the location 
of the first peak well, though the spacings between other 
peaks is better. The above values were calculated assuming 
h = 0.65, Sib = 0.05, n = 1 and ai^ = 1100"^ 



which we measure today. The odd peaks correspond to 
maximum compression of the fiuid, the even ones to rar- 
efaction [10]. In an idealised model of the fluid, there 
is an analytic relation for the location of the m-th peak: 
l^ ^ uiIa [11, 12] where Ia is the acoustic scale which may 
be calculated analytically [6] and depends on both pre- and 
post-recombination physics as well as the geometry of the 
universe. 

The simple relation Im ^ ttiIa however does not hold 
very well for the flrst peak (see Table 1) although it is bet- 
ter for higher peaks [2]. Driving effects from the decay of 
the gravitational potential as well as contributions from the 
Doppler shift of the oscillating fluid introduce a shift in the 
spectrum. In order to compensate for this, we parametrize 
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the location of the peaks and troughs as in [11] by^ 

Im =lA{m- ifm) =lA{m- (p - 5<fim) ■ (1) 

For convenience, we define (fi = ipi to be the overall peak 
shift, and 5ipm = — (p the relative shift of the m-th peak 
relative to the first. The reason for this parametrization is 
that the phase shifts of the peaks are determined predom- 
inantly by pre-recombination physics, and are independent 
of the geometry of the Universe. In particular, the ratio of 
the locations of the first and m-th peaks 

\rn_ _ Ia {m- (fi - S'^rn) _ ^ TW - 1 - Sifm 

h ~ Ia (1 - ^) ~ ^ 1 - ^ ' ^ ' 
probes mostly pre-recombination physics and so can be used 
to extract information on the amount of dark energy present 
before last scattering [6]. 

If wc knew how the phase shifts depended on cosmologi- 
cal parameters, it would be possible to extract Ia from the 
measured CMB spectrum. Since any given cosmological 
model predicts a certain value of I a, this is a simple way 
of distinguishing between difi^erent models - in particular it 
has been shown [6] that different quintessence models with 
the same energy density and equation of state today can 
have significantly different values of I a- Finally, having ex- 
tracted Ia from observations, we could speed up likelihood 
analysis by being able to discard models not leading to the 
right value of the acoustic scale before a single perturbation 
equation has to be solved. 

In a recent paper [11], a fitting formula for was given 

^« 0.267 (^) , (3) 

for the values n = 1, Qbh^ = 0.02. In this formula, r* is 
the ratio of radiation to matter at last scattering^ 

r* = Pr{z*)/Pm{z*) = 0.042 (Qmh') (0*/lO^) . (4) 

Equation (3) however, is valid only for the given values of 
spectral index, Hubble parameter and baryon density. It 
does not include the dependence of the peak location on 
the amount of quintessence present at last scattering, and 
is valid only for the first peak h. In this paper, wc give 
fitting formulae (sec Appendix A) for the shifts of the first 
three peaks and the first trough and describe how one can 
use them to extract cosmological information from future 
CMB experiments. 

Our first task in computing fitting formulae for the peak 
locations is to decide which cosmological parameters to fit 
to. The dependence on the baryon density and the Hub- 
ble parameter is sensitive only to the product Qbh!^, and 
so we do not seek to fit for them separately. We further 
take r* defined in Equation (4) and the spectral index n as 
parameters. For the quintessence dependence, we use the 
effective average density component before last scattering 
fiit defined as in [6] 

n,t = T,-' r n^r)dr. (5) 

^0 

^The peaks are labelled by integer values of m and the troughs 
by half-integer values. 

^This relation also holds in the presence of quintessence. 



Symbol 


Range 




[0.2, 0.6] 




[0.005, 0.04] 


Hit 


[0, 0.23] 


h 


[0.55, 0.80] 


n 


[0.8, 1.2] 



Table 2: Parameter ranges used in this work. 



We recall that the peak shifts are sensitive mainly to 
pre-recombination physics and so we do not need to use 
the value of fl'* today a.s a parameter. Of course the acous- 
tic scale Ia does depend on today's quintessence component 
- we give a relation for Ia in Section 3. We will thus seek 
to find the dependence of (<^, 5ip,n) on the cosmological pa- 
rameter set (^flih^ ,r* ,n,flfs^ . In performing these calcu- 
lations, we restricted each of the cosmological parameters 
used to lie within a certain interval, which in each case is 
over- rather than under-cautious. The ranges of parameter 
values chosen are displayed in Table 2. To gain intuition 
for the fitting formulae, we plot curves for the shift of the 
first and the second peak as well as the relative shifts of the 
first trough and the second peak in Figure 1. 

In Sections 2 and 3 we describe a systematic procedure 
for extracting the acoustic scale I a from the location of the 
first three peaks. Section 4 introduces a quantity n which 
is useful as it depends only on two of our four parameters. 
The model (in) dependence of the fitting formulae is dis- 
cussed in Section 5. Finally, our fitting formulae are given 
in Appendix A. 

2 Retrieving the shifts from 
CMB measurements 

With future high precision measurements of the MAP^ and 
PLANCK^ satellites, we expect that the position of the first 
three peaks and troughs will be determined to high accu- 
racy. From these few data points, it is possible to extract 
valuable information on the cosmological parameters. We 
have observed, during our computation of CMB spectra for 
thousands of universes, that the overall shift of the third 
peak (p3 (i.e. ips = ip + Sip-^,) is a relatively insensitive quan- 
tity. In the parameter range wc used (see Table 2) we found 
that (p3 = 0.341 ± 0.024."'' In using (p3 — 0.341 we introduce 
slight (at most one percent) systematic deviations in our 
estimate, because an increase of f2is typically increases ips 
(see Fig. 1(d)). We will partially correct for these eff'ccts by 
improving our estimate for (ps , via the procedure described 
below. 



^http://map.gsfc.nasa.gov/ 

*http://astro.estec.esa.nl/SA-general/Projects/Planck/ 
^Here and in the following, we quote 1-a errors. All errors 
follow approximately a bell curve. 
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Figure 1: The overall shift i^s (a) and the relative shifts of the first trough (b) and the second peak (c). Also given 
is the overall shift of the third peak (d). In all figures, the long dashed, dotted and the dashed lines represent 
the fitting formulae for the parameters ^f2f,/i^, r*, n, iljt^ = (0.02, r*, 1, 0), (0.02, r*, 1, 0.1) and (0.01, r*, 1, 0) 
respectively. The large symbols show the data corresponding to these curves. The errors quoted in Appendix A 
are calculated from the spread of these symbols relative to the curves. The sprinkled dots represent thousands 
of models selected at random from the parameter space given in Table 2, and indicate the ranges of values taken 
on by (p etc. for these models. 



3 





(%) 


/, _num\ 


^^improvod ^ 


0- 


2 


0.313 


0.326 


10- 


12 


0.340 


0.337 


18- 


20 


0.362 


0.348 



Table 3: Binned average ips of the numerical simulation and 
the improved deduction. 



We start by extracting our first estimate of the overall 
phase shift, from the measured locations of the first and 
third peaks 



(3- 



Ji 



2.66f . 

<3 



(6) 



Comparing this estimate with the value calculated from nu- 
merical simulations, we find Aip — 0.006. Having a handle 
on the overall phase shift, it is now simple to infer the rel- 
ative shifts 5ipm of the remaining troughs and peaks. From 
equation (2) we get the relation 



Sifim = (m - 1) 



The error of this estimate is 



(7) 



(8) 



Having a first (and already quite accurate) estimate of the 
shifts, we now correct for the systematic effects described 
above. Taking the cosmological parameter set we wish to 
maximise over (i.e. Table 2), we calculate for each model 
universe the phase shifts of the first three peaks using the 
fitting formulae given in Appendix A. We then discard 
those models for which any phase shift deviates significantly 
(say > 2-a) from the data-inferred values. This leaves an 
improved cosmological parameter set, for which the aver- 
age value of (p3 is calculated (see Table 3). This improved 
ifis can then be used to re-calculate the phase shifts from 
Equations (6) and (7). 



3 Estimating Ia 



Using the improved value® for (ps from the previous section, 
we can extract to very good accuracy the acoustic scale I a - 
the quantity which determines the overall spacing of CMB 
peaks: 



Ia = 



h 



(9) 



In fact, the deviation of the value of Ia estimated from this 
formula and the rmmerically-obtained value is small, with a 
1-cr error of 0.8% (see also Table 4). This is a very valuable 
result, for the value of I a can be simply computed for any 



®In fact, using <p3 = 0.34 instead of the improved value also 
gives reasonable results. 



given quintessence (or indeed any other) cosmology. For 
flat universes it is given by [6] 



Ia 



— 1 

TTC, 



F(f]^,Wo) 



Ols -I- 



1-n^ 



(10) 



with 



3Wo) 



05(1 -g) 



-1/2 



(11) 



where D,q,D,q are today's radiation and quintessence com- 
ponents, ais is the scale factor at last scattering (if ao = 1) 
and Cs is the average sound speed before last scattering: 

c. = Til^ dr [3 + i9/4){p\t)/p-'{t))] . (12) 

The effective equation of state. Wo is the fi'^-weighted av- 
erage over conformal time 



Wo = 



n^(r)«;(r)dr X fi'^(r)drj . (13) 

In particular, different quintessence models with the same 
energy density and equation of state today can have signif- 
icantly different values of I a- In this way stringent bounds 
on cosmological models can be imposed just by comparing 
the Ia value of specific models. 



4 Insensitive Quantities 

The phase shifts depend on the cosmological parameters 
course, if it were possible to find a 
linear combination of phase shifts which is insensitive to 
some of these parameters and thus reduce the dimension- 
ality of our parameter space, it would greatly help in ex- 
tracting cosmological information. To this end, we note an 
anti-correlation between ip and Sips - empirically, we have 
found that the quantity 

2 

K=ip+ -Slfi3 (14) 
o 

is practically insensitive to and flth^, and depends only 
on n and fii^. In fact, it is to very good approximation 
given by the fit 

K = (0.277 + 0.284fiit) (1.3 - 0.3n), (15) 

with Ak*** Ri 0.0024 being the deviation of the fit from the 
immcrically-sinmlatcd values (see Fig. 2). Following the 
procedure in Section 2, we can deduce k from the measured 
values of the peak locations. Within our parameter range, 
K is then determined with error Ak'^'"^""- = 0.013. 
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0.2C - 



0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 

Figure 2: The quantity function of r*. It is practically 

insensitive to and ilt/i^ for most of the initial conditions 
considered. The dots represent fifty thousand models with 
parameters in the ranges given in Table 2 The +'s and x's 
represent models with D.fs = and 0.22 respectively, for 
n = 1, and all values of other input parameters. 



In the parameter space we have considered, the value of 
K varies between 0.26 and 0.36. Hence to 1-a confidence 
level, about three quarters of our two-dimensional (n, fiit) 
parameter space can be excluded for any given k. For in- 
stance, without quintessence, the value of k lies between 
0.26 and 0.29 for n € [0.8, 1.2]. The measurement by MAP 
or PLANCK of a value of k > 0.29 would therefore be a 
strong hint of a dark energy component playing a role at 
last scattering. 



5 Model dependence 

The fitting formulae were obtained using a standard expo- 
nential potential [13] for the quintessence component. Be- 
cause the shifts are almost independent of post recombi- 
nation physics, we expect the results to be approximately 
correct for any realization of quintessence, i.e. all poten- 
tials. One should however be cautious with models that 
are qualitatively extremely different from the exponential 
potential before last scattering, as for example the Ratra- 
Peebles inverse power law [14] with substantial ilf^. In these 
models there is a sharp increase in Q,'^ during recombina- 
tion, whereas the quintessence content for the exponential 
potential is fairly constant at this epoch. 

The inverse power law is characterized by its potential 
yiPL = A/ip'^. Models with a ^ 2 are phenomenologically 
disfavoured [15]. We use these models only as cross checks 
for the fitting formulae. 

In terms of phase shifts, one finds that the sensitive rel- 
ative shifts of the first trough and the second peak diff'er 
substantially for the two models (see Table 4). However, (f 
and K are seen to be more robust and the deduced value of 
Ia is accurate to within one percent in every case. 



6 Conclusions 

In this paper we have shown that within a wide range of 
parameters, one can accurately deduce the acoustic scale 
I A, as well as the shifts of the peaks and troughs provided 
the locations of the first three peaks are measured. Not 
only will this enable faster testing in likelihood analysis 
by providing a filter before any fiuctuation equations are 
solved, but it could also in principle lead to a detection 
of quintessence - measuring a non-zero value of the dark 
energy at last scattering (e.g by computing the quantity k 
described in Section 4) would distinguish it from a cosmo- 
logical constant, whose contribution to the energy density 
of the Universe would become significant only very recently. 
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A Fitting formulae 

We present here our fitting formulae for the overall phase 
shift if, followed by the relative shifts of the first trough 
(3^3/2) and the second {S1P2) and third (Sips) peaks. ^ In 
each case we also give an estimate of the accuracy of the 
formulae. 

A.l Overall phase shift Lp 

For the overall phase shift (f (i.e. the phase shift of the first 
peak) we find the formula 

>f = (1.466 - 0.466n) [air^^ + 0.29init] , (16) 
where ai and a2 are given by 

ai = 0.286 0.626 (STib/i^) (17) 

a2 = 0.1786 - 6.308 06/1^ + 174.9 (S^b/i^)^ (18) 

-n68{nth'f- (19) 

It contains the main dependence of any shift iprn on D,f^. 
The 1-a error for ip is 

A.^ = 0.0031 (20) 

A. 2 Relative shift of first trough 6(^^/2 

The relative shift of the first trough is a very sensitive quan- 
tity spanning a wide range of values. It can very well be 
used to restrict the allowed parameter space for cosmolog- 
ical models. We have 

5^3/2 = feo + bir-y^ exp(62r*) + 0.158 (n - 1), (21) 

' A small cH — h package providing functions for the 
shifts is available at http://www.thphys.um-heidelberg.de/~ 
doran/peak.html 
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Ols (%) h I3/2 h h Ia V 3(p3/2 5ip2 K 



Leaping kinetic term 



3 


214 


396 


521 


788 


293 


0.269 


-0.121 


-0.045 


0.287 












294 


0.271 


-0.119 


-0.041 


0.292 


13 


210 


396 


522 


799 


301 


0.301 


-0.120 


-0.038 


0.317 












301 


0.301 


-0.120 


-0.038 


0.318 


22 


208 


397 


524 


808 


307 


0.324 


-0.116 


-0.030 


0.341 












305 


0.320 


-0.120 


-0.035 


0.333 






Ratra Peebles 


inverse power 


law 






: 10"^ 


199 


366 


480 


724 


269 


0.259 


-0.119 


-0.043 


0.278 












270 


0.261 


-0.117 


-0.038 


0.284 


10 


178 


339 


443 


674 


251 


0.294 


-0.140 


-0.054 


0.304 












253 


0.298 


-0.138 


-0.050 


0.312 


22 


172 


338 


444 


683 


258 


0.333 


-0.144 


-0.057 


0.340 












258 


0.334 


-0.145 


-0.057 


0.340 



Table 4: The peak locations and the phase shifts of leaping kinetic term [16] and Ratra Peebles inverse power law [14] models 
for Qbh'^ = 0.021, Q,^ = 0.6, h = 0.65, n = 1 and varying file. The inverse power law models correspond to a = 6, 22 and 
40 respectivly. The first row of each model gives the CMBFAST-obtained values of the locations of the peaks and the phase 
shifts as well as Ia and k. The second row gives the values deduced using the method described in Section 2. 



with 



60 = -0.086 - 0.079 nt ~ (2.22 - IS.mt) Q.bh^ 



A. 4 Relative shift of third peak 5ifz 

For the third peak, we find 





- [l40 403f7it] {Q.bh^ f 


(22) 


6(fi3 = 10 — dirf^ 


+ 0.08 (n 


61 


= 0.39 - 0.98nit - (18.I - 29.2ait) fifc/i^ 


(23) 


with 






+440 (Qbh^)^ 


(24) 


di = 9.97 -1- (3.3 - 3Ql' 


) Qbh'^ 


62 


= -0.57 - 3.8 exp |-2365.0 (fifc/i^)^! . 


(25) 


d2 = 0.0016 -0.0067 af^ 


, + (0.196 


For the 


one standard-deviation error we have 




(2.25 2.77 nfj 
and error given by 


X 10"^ 




ASip3/2 = 0.0039. 


(26) 




A.3 


Relative shift of second peak Sip2 


A6ip3 = 


0.0052. 



The relative shift of the second peak is a very sensitive 
quantity. It is thus not surprising to find a strong depen- 
dence of 6(p2 on the parameters. We have 

5<fi2 = co- cm - C2r-=^ + 0.05 (n - 1), (27) 

with 

Co = -0.1 (^0.213 - 0.123 Hit) (28) 

X exp |- (52 - 63.6 Hit) fib/i^ j (29) 

ci = 0.063 exp I -3500 (J^b/i^)^ I 0.015 (30) 

C2 = 6 X 10"" 0.137(^6/1^ -0.07)^ (31) 

C3 = 0.8 + 2.3 nt + (70 - 126 Hit) flbh^- (32) 
The error of this approximation is 

A6ip2 = 0.0044. (33) 



(34) 

(35) 

(36) 
(37) 



A. 5 Overall shift of third peak cp^ 

For completeness, wc give a fit for ip3 which in principle 
could be obtained by adding (p and Sips. However, a one- 
step-fit yields better errors here. Our formula is 



953 = ei (1 -h e3r*) rl^ + a - 0.037 (n - 1), (38) 



with 



ei = 0.302- 2.112 na/i^-h 0.15 exp {-384^6/1^1(39) 
-0.04-4.5^6/1^ (40) 

(41) 



62 



63 = (-0.118 + 44.7 06/i^) fiit 

64 = (0.214exp {-48^6/1^} +0.106) fiit, (42) 
and error 

A(p3 = 0.0017. (43) 
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